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ABSTRACT 

We present near-infrared J, H , and X images of four embedded stellar clusters in the Galaxy. 
We find a significant fraction of pre-main-sequence stars present in at least one of the clusters. For 
the clusters dominated by main-sequence stars, we determine the initial mass function (IMF) both 
by using the K luminosity function and a global extinction correction and by deriving individual 
extinction corrections for each star based on their placement in the A" vs. H — K color-magnitude 
diagram. Based on our IMFs we find a significant discrepancy between the mean IMF derived via the 
different methods, suggesting that taking individual extinctions into account is necessary to correctly 
derive the IMF for an embedded cluster. 

Subject headings: open clusters and associations: general — stars: formation — stars: luminosity 
function, mass function 



1. INTRODUCTION 

Embedded clusters are increasingly recognized as vi- 
tal sites of star formation for both low- and high-mass 
stars. Recent studies indicate that clusters may ac- 
count for 70-90% of star formation and that embed- 
ded clusters (those still partially or fully enshrouded 
in their natal molecular cloud) may exceed the num- 
ber of olde r, non-embedded open clusters by a fac - 
tor of -20 <|Elmegreen et al.ll2000t lLada fc LadH kOO.T). 
The stellar content of embedded clusters within well- 
known star formation regions can now be probed, where 
high extinction (Ay > 10) prohibits studies at opti- 
cal wavelengths. The IMF of such clusters has gener- 
ally been found to be consistent with a Salpeter value 
with a slope of T = —1 - 35 (e . g lOkumur a et alJ 1200(1 
iBlum. Damineli. fc ContH I200H iFigueredo et alJ 12002ft 
although outliers have been found as well, g enerally with 
flatte r slopes than the Salpeter value (e.g. iPorras et al.1 

Although near-infrared (NIR) spectral classi fication of 
massive stars is possible ijHanson et al.llT996ft . in most 
cases determinations of the IMF from NIR data rely 
heavily, if not exclusively, on photometry and use spec- 
troscopy only to obtain reliable masses of the few bright- 
est and most massive stars in a cluster if at all. Since 
these results thus depend on stellar evolutionary models 
as well as details of the handling of extinction, this raises 
concerns about to what extent the IM F depends on the 
methodology employed. [Massevi $2003j) cites the example 
of N GC 6611, where two separate analyses of the same 
data ({Hillenbrand et aT]ll993i iMassev et al.lll99"5ft using 
different treatments of extinction produced IMFs differ- 
ing by more than the formal la errors would suggest 
(r = -1.1 ±0.1 and T = -0.7 ±0.2.) Similarly, the IMF 
for the G305±0.2 embedded cluster differs by more than 
the errors between the value derived from the K lumi- 
nosity function (KLF; T = —1.5 ± 0.3) and that derived 
using the color-magnitude diagram (r = —0.98 ± 0.2) 
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ijLeistra et al.l l2005). Claims that variations in the IMF 
exist, whether based on individual extreme cluster s such 
as th e Arches or a general analysis of the data ijScald 
1998), must thus be handled with care to compare only 
results based on similar methodology. 

The final release of the Two Micron All Sky Sur- 
vey (2MAS S) has fost ered s tudies (e.g. IDutra fc Bical 
l2000l l200lt iBica et alJl2003 llvanov et al] 12002ft which 
can probe a much larger portion of the Galaxy for pre- 
viously unknown embedded stellar clusters and signifi- 
cantly increased the number of known embedded clus- 
ters. A compilation of some of these results along 
wi th previously known embedded clusters is presented 
in IPorras etldl l|2003fl who find that - 80% of the stars 
in their sample are found in "large clusters" of more than 
100 stars, despite the rarity of such clusters. However, 
these studies are not foolproof, and compilations of "em- 
bedded clusters" based purely on the 2MASS data with- 
out followup must be treated with cautio n. The studies 
based solely on stellar density criteria (e.g. IDutra fc Bical 
2000, 2001) have been found in fol l owup work by dif- 
ferent groups (e.g. IDutra et al.ll2003t iLeistra et al.ll200"5l 
iBorissova et al.ll2005ft to have only about a 50% success 
rate toward the inner Galaxy where the stellar back- 
ground is high. We have performed an independent 
search of the 2MASS archive, searching the Point Source 
Catalog for regions of higher stellar density than the 
background (determined locally within a 5' radius) which 
are redder in H — K than the local field. This selects 
for embedded clusters, with the color criteria helping to 
eliminate chance superpositions and regions of low ex- 
tinction. A large background radius and the use of color 
selection are critical to the automated identification of 
embedded clusters, but even color selection can fail in 
regions of high background stellar density, where patchy 
extinction can mimic clusters. 

In ||21we present the observations and data reduction 
for four embedded clusters found in the 2MASS Point 
Source Catalog, in &\ we present the A-band luminos- 
ity functions (KLF) and initial mass functions (IMF) for 
the clusters, and in H3.5l we address the issue of system- 
atic differences between different methods of deriving the 
lalMfeonS'edipkedded clusters, for our clusters as well as 
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IMFs from the literature. 

2. OBSERVATIONS & DATA REDUCTION 

We selected five young stellar cluster candidates from 
the 2MASS Point Source Cat alog based on color and 
density criteria as described in ijLeistra et aUl2005|) . We 
selected regions with a higher stellar density than the 
locally defined field with redder H — K color than the 
field to select embedded cluster candidates. The cluster 
candidates were observ ed using the PISCES instrument 
ijMcCarthv et al.ll2001j) on the 6.5m MMT on JanlO-11, 
2003. PISCES uses a 1024x1024 HAWAII array with a 
platescale of O'.'185/pixel on the MMT, providing a 3'x3' 
field of view. Images of all cluster candidates were ob- 
tained in J, H, and K filters to a limiting magnitude of 
J = 19.5 , H = 18.5, K = 18. Four of our candidates (all 
except the one near Sh 2-21 7) were independe ntly identi- 
fied as cluster candidates bv lBica et alJ l)200"3j) . who used 
criteria based only on stellar density without color con- 
siderations. Our deeper, higher-resolution images sug- 
gest that four of these candidates are genuine stellar 
clusters, while the fifth, near Sh 2-258, contains only a 
few stars in the PISCES images and could be either a 
very small cluster or a chance superposition. We present 
our results for the four confirmed clusters in this paper. 
Seeing conditions when the images were obtained were 
variable, ranging from 0'.'5 to l'.'l at K . 

All images were reduced and combined using IRAF 
routines. The distortions of the PISCES camera were 
mapped by imaging the globular cluster NGC 4147 and 
mapping the observed locations to the USNO-B known 
coordinates, then constructing a transformation function 
using the IRAF task geomap and correcting the images 
using IRAF geotran. There are too few USNO-B stars 
in the heavily extincted regions we observed to provide 
suitable distortion corrections from the fields themselves. 
The individual images, taken with a spiral dither pattern, 
were then combined. We fit a PSF to each image using 
the IRAF task psf, allowing the PSF to vary across the 
field to compensate for residual distortions. Photometry 
was then done using IRAF-DAOPHOT. 

Photometric calibration was performed using the 
2MASS magnitudes of field stars. Stars used for calibra- 
tion were selected to be well-separated from other stars 
and from nebulosity in the PISCES images to ensure 
that they were uncontaminated in the lower-resolution 
2MASS images, and to have magnitudes bright enough 
to have good photometry in 2MASS (K < 14). We 
chose to use a relatively large number of calibration stars 
rather than selecting the few most isolated stars to reduce 
effects of potential variability and photometric outliers 
among the calibration stars. The scatter in the photo- 
metric calibration derived from comparison to 2MASS is 
the dominant source of photometric error, contributing 
two to three times the measurement errors as reported 
by DAOPHOT. DAOPHOT errors were ~ 0.03 magni- 
tudes while the calibration uncertainties were ~ 0.1 mag- 
nitudes. Quoted errors in the 2MASS photometry were 
negligible, with most stars having an error of ±0.003 mag 
or less in all bands. Thus, the quoted error should be 
considered an overestimate when considering the relative 
photometry of stars within either cluster; the calibra- 
tion errors from comparison to the 2MASS photometry 
will shift all our measurements by the same amount. No 



trend with location on the chip was observed in the cal- 
ibration for any of the clusters, though the scatter be- 
tween the PISCES and 2MASS magnitudes becomes sig- 
nificant in the outermost 15 "; we thus exclude these 
sources from the analysis. 

3. ANALYSIS 

3.1. Sh 2-217 Cluster 

We present a if-band image of the cluster near Sh2-217 
in Figure n The cluster is nearly circular in projection 
and is quite dense; even in the highest-resolution indi- 
vidual pointings we obtained it suffers from crowding in 
the central regions. The cluster extends over most, if not 
all, of the field of view. We present the K image rather 
than a color frame because the seeing was significantly 
better at K. 

This cluster was analyzed in the NIR by 
iDeharveng et alJ l)2003f) . who discuss the large un- 
certainties in the distance to Sh2-217. Based on the 
Lyman continuum fluxes from the main exciting star of 
Sh2-217 (located several arcminutes outside our field of 
view; the cluster is located on the periphery of the H II 
region) they adopt a distance of 5.0 ± 0.8 kpc, which is 
consistent with the kinematic distance to the associated 
molecular gas. The cluster is coincident with a peak in 
the 8 jum emission as measured by the MSX mission, 
suggesting dust is still present in the cluster. 

The K vs H — K color-magnitude diagram is shown 
in Figure El The stellar density (determined in K where 
the seeing was best) does not plateau in the 3' x 3' field of 
view, suggesting that the true field-star density level has 
not been reached, and cluster stars are still present out to 
the edges of the field. As a result, in order to correct for 
foreground contamination, we selected an adjacent field 
of the same size from 2MASS to use as a comparison 
field. We assumed the luminosity function of the field to 
be the same as that of the outer portions of the cluster 
(excluding the inner regions to minimize the results of 
crowding and mass segregation) in order to extrapolate 
from the limiting magnitude of 2MASS to that of our 
images. We then binned the field stars by K and H — K 
with a bin size of 0.5 magnitudes and randomly selected 
the appropriate number of stars for removal from each 
bin in the cluster region. Thi s is similar to the p rocedure 
employed by, among others, (B lum et al.ll2000|) (who do 
not describe an extinction correction; this lack of correc- 
tion is equivalen t to as suming a common extinction) and 
iFigueredo et al J 1 2002 ). and is the method we employed 
in l)LeisteTet*liLl 2005). The resulting statistically cor- 
rected CMD is shown in Figure A total of 62 stars 
were removed in this procedure, out of an initial total of 
236. The fairly wide distribution in H — K for cluster 
stars is likely due to a combination of factors, notably 
an actual spread due to differential extinction to differ- 
ent regions of the cluster and to the greater influence of 
crowding in H (where the seeing was poorer). Although 
we find that individually correcting extinctions gener- 
ally provides a superior estimation of the IMF compared 
with using only the K data and a single extinction for 
the cluster as a whole, in this situation the lower quality 
(in particular the poorer seeing and consequently more 
severe crowding) of the 7J-band data leads us expect that 
the K luminosity function (KLF) will produce a better 
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estimate of the IMF for this cluster than the CMD will. 
We have previously compared these two m ethods of de- 
termin ing the IMF for embedded clusters in lLeistra et al.l 
(2005); in that case, we found they gave different results, 
with the "CMD" method (which we anticipate will be 
more reliable in most cases, especially where variable ex- 
tinction is present) yielding a flatter slope. 

3.1.1. The KLF 

In order to obtain a robust determination of the KLF 
and the IMF, we need to determine the completeness of 
our data. To do this, we performed artificial star tests. 
We inserted five artificial stars, each of the same magni- 
tude, at a time into the cluster region, then ran IRAF- 
DAOPHOT with the same parameters as we used for the 
initial analysis. This procedure was repeated 50 times for 
each magnitude bin (Am = 0.5), for a total of 250 arti- 
ficial stars added in each bin in H and in K. The stars 
were added in small numbers at a time to avoid hav- 
ing the artificial stars significantly change the crowding 
characteristics and thus influence the completeness. The 
artificial star tests indicate a high level of completeness 
down to K = 17.5. The actual completeness is most 
likely slightly lower, since in the crowded central region 
of the cluster our method may produce false positives 
when the artificial star is placed on top of a real star of 
approximately the same magnitude. Despite this concern 
we have used the calculated incompleteness in correcting 
the KLF; however, we have excluded the K = 17.5 bin 
from consideration, both because this effect will be most 
pronounced at faint magnitudes, and because statistical 
uncertainties in the incompleteness will be significant. 

Knowing our incompleteness, we can calculate the KLF 
for the cluster. Figure^shows both the uncorrected and 
completeness-corrected versions of the KLF, as derived 
from all sources detected in K . The slope of the KLF is 
0.31 ± 0.04, with no extinction correction applied. 

3.1.2. The IMF 

We derived an IMF for the Sh 2-217 cluster by two 
methods. For bot h methods we used a d istance to the 
cluster of 5 kpc ijDeharveng et alJ l2003|) . In general 
we believe that the "CMD method" for deriving the 
IMF, which uses individually-derived extinctions for each 
star in the cluster, to be more reliable than the "KLF 
method" which assumes a common extinction to all stars 
in the cluster, since variable extinction is frequently ap- 
parent in the NIR images of embedded clusters. How- 
ever, in this case our K data is superior to our H data 
due to the difference in the seeing, which was ~ 0." in K 
and ~ 1.2" in H . This suggests that despite the general 
drawbacks of the KLF method it may be preferable in 
this situation; including the -ff-band data adds nothing 
if it is of poor quality. At the least, using both methods 
will provide information on potential systematic effects 
in the IMF determination that depend on the methodol- 
ogy used. 

The KLF method of determining the IMF is some- 
times employed even whe n multi-color pho tometry and 
spectra are available (e.g. iBlum et al.ll2000|) . and is sim- 
ply a transformation from K magnitude bins to mass 
bins. The major problem with this is that of variable ex- 
tinction, which for many embedded clusters is significant 



even in K . To make this transformation, we first correct 
the observed K for distance and extinction. Without 
spectra, we cannot obtain a precise estimate for the ex- 
tinction; instead, we compute an average extinction cor- 
rection based on the observed J — H and H — K colors of 
the brighter stars. Since there is little difference in the 
intrinsic H — K color of stars of spectral type F5 and 
earlier, this estimate of an average extinction is not sen- 
sitive to minor errors in the di stance estimate. Us i ng the 
stellar evolutionary models of iMevnet &: Maederl (2003) 
for solar metallicity, we relate the mass for each track 
to an absolute K magnitude for a star on the ZAMS. 
We trans formed Lhni to K using the bolom etric correc- 
tions from lVacca. Garm anv. fc Shul| (119961) for the early 
spectral types anc HMldligninTet*ld.l l)1986j) for later spec- 
tral types. We then us e the intrinsic V — K col ors from 
iBessell fc Brett! l)1988j) for A-M stars and from iWegnerl 
(1994) for O and B stars. Finally we interpolate linearly 
between the masses available on the evolutionary tracks 
to find the masses corresponding to our magnitude bins, 
and fit a power law to the resulting mass function. The 
IMF slope we derive by this method is T = -2.7 ± 0.25, 
excluding bins corresponding to K > 17.5 where in- 
completeness becomes significant. The slope we fit to 
the KLF itself for sources detected in both H and K is 
0.35 ±0.04. 

We have previously used multi-color photometry in 
conjunction with near-IR spectroscopy of the bright- 
est stars to determ ine the IMF for embedded clusters 
<|Leistra et all 12001 . Although we do not have spec- 
tra in this case, we can still derive extinctions for in- 
dividual stars based on their near-IR colors. We use 
the same evolutionary models and conversions from the- 
oretical to observed quantities described for the KLF 
m ethod. The ZAMS d e rived from the evolutionary tracks 
of IMevnet fc Maederl {2003) is ovcrplotted on the dis- 
tance and extinction-corrected CMD in Figure 03 The 
ZAMS lies in the middle of the distribution of stars due to 
the method used to estimate an average extinction. The 
scatter around the ZAMS is rather large, and is likely 
due to a combination of variable extinction in the clus- 
ter region and poor photometry in H, especially in the 
central portion of the cluster. 

Since extinction appears to vary across the cluster, we 
impose an extinction limit on the sample used for the 
CMD computation of the IMF. Since the most massive 
stars can be seen to greater extinction than less massive 
stars, neglecting to impose this constraint will produce 
an overly flat IMF. Thus we need to simultaneously limit 
our sample by extinction and by mass. We use a mean 
extinction to the Sh2-217 cluster of Ay — 9 mag as de- 
termined by individually de-reddening sources until they 
reach the main sequence. With the exception of a few ex- 
treme outliers that most likely suffer from poor photom- 
etry, most sources with higher extinctions have Ay < 15. 
At a distance of 5 kpc with our limiting magnitudes, we 
can observe stars earlier than G4 to an extinction of 9 
mag and earlier than F0 to an extinction of 15 mag. 

At a distance of 5 kpc with our limiting magnitude, an 
extinction of Ay — 9 limits us to G4 and earlier stars, 
while Ay — 15 limits us to F0 and earlier. We select 
Ay — 10 and G2 as our limits; stars with higher ex- 
tinction or later spectral type cannot be detected over 
the entire range of mass or extinction included and thus 
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are excluded. Approximately 34 stars have a higher ex- 
tinction than this, including five with extreme calculated 
extinctions (Ay > 50) that most likely suffer from poor 
photometry and have unrealistic colors. This extinction 
limit will also have the effect of excluding stars with 
if-band excess from the IMF determination, since such 
sources would appear to be at high extinction. This will 
tend to push the IMF to flatter values, since lower-mass 
sources spend more time as IR-excess objects and thus 
are more likely to be ruled out by this criterion. However, 
we consider this to be a better approach than including 
the sources since: 1. the number of sources detected in 
all three bands in this cluster showing near-IR excess is 
small, suggesting such sources will not significantly in- 
fluence the mass; 2. the majority are of a low enough 
mass to fall below the completeness limit, and are thus 
excluded anyway regardless of the method; 3. includ- 
ing them (since not all sources are detected in J) would 
weaken the extinction limit and tend to again force the 
IMF to flatter values. In clusters where IR-excess sources 
dominate we do not fit an IMF (see Section 13. 2. II for fur- 
ther discussion). 

Individual extinctions are derived by moving the stars 
along the direction of the reddening vector until they lie 
on the ZAMS. Once the stars have been corrected indi- 
vidually for extinction, they are placed in mass bins. We 
fit a power law to the data, excluding masses < 1.1M Q 
from the fit since they cannot be seen over the entire 
range of extinctions in the cluster. The IMF slope we 
derive by this method is T = —1.61 ± 0.2. As for the 
KLF method, the quoted errors represent only the for- 
mal errors in the fit and should be considered an under- 
estimate. 

The difference between these two values for the IMF 
emphasizes that formal statistical errors significantly un- 
derestimate the true un certainties in the IM F. In this 
case, as was the case in ILeistra et all ll200l . the CMD 
method gives a noticeably flatter result than the KLF 
method. Clearly the individual extinction correction 
leads to the conclusion that more massive stars are 
present than an average correction does. This could 
be due to the effects of mass segregation, or to an in- 
correctly chosen extinction limit (so that we truly are 
sampling massive stars more completely than lower mass 
stars). It is difficult to understand which of these ef- 
fects is most important without obtaining spectra for a 
significant number of stars in the c luster. 

This cluster is also analyzed bv iPorras et alJ 0999), 
who use a slightly different distance (5.8 kpc) and ex- 
tinction (< Ay >= 5.3 ±3.7) and derive an IMF slope of 
T = —0.59 based on 54 sources using the J vs J—H CMD 
to individually correct extinctions and compare with a 
theoretical JLF. This is a significant discrepancy from 
our result with either method. A number of factors may 
contribute to this difference. The most significant, how- 
ever, is likely to be due to a different choice of cluster 
boundaries. Their quoted cluster radius corresponds to 
only 50'.'9, compared to ours of ~ 80". This suggests 
that their IMF will be more weighted toward the cluster 
core than ours. The value they quote for a field + cluster 
IMF is r = —2.71 ± 0.24, much steeper and in fact quite 
close to the value we obtain using the KLF. If the cluster 
suffers from mass segregation, as we would expect given 
that it is observed even in very young clusters (e.g. the 



Arches cluster (jStolte et al.ll2002|) L we expect the core 
IMF to be quite flat. When we re-derive the IMF for 
this cluster using their radius with our data, we derive 
an IMF of T = —1.55 ± 0.22 , statistically indi stinguish- 
able from our original result. IPorras et alJ l)1999fl do not 
comment on issues of confusion or field star contamina- 
tion, so we cannot evaluate how much of an effect it is 
likely to have on their result; we expect crowding to be a 
more significant issue, and a misidentifi c ation of blended 
sources as single stars bv IPorras et alJ 0999) could ac- 
count for their finding a steeper IMF than we do using 
the same method. 

3.2. IRAS 06058+2158 Cluster 

We present a three-color composite of the cl uster near 
the IR AS source 06058+2138 in Figure El iBik et alJ 
( 2005) obtained VLT spectra in K of several NIR point 
sources near the IRAS point source, which is located 
near the center of the cluster, and identified two can- 
didate massive YSOs and an embedded early-B star. 
The spectrophotomctric distance they derive from the 
B star is 1.0-1.5 kpc. This cluster is much more heav- 
ily embedded than the Sh 2-217 cluster, with significant 
nebular emission and prominent dust lanes. Numerous 
OH and methanol masers have been detected in this re- 
gion, which along with the IRAS point source sugges t 
ongoing star formation , (see, e.g.. iCaswell et all (1995), 
iSzvmczak et alJ {2Q00)). A peak in the 8 yum emission 
is observed in the MSX data, extending from the region 
of NIR nebulosity to the southeast to the isolated bright 
source. 

The embe dded cluster h e re is described in the com- 
pilation of lLada fc Ladal il2003l). who q uote a dis- 
tance of 1.5 kpc llCarnenter et all 119931) . However, 
lHanson. Luhman. fc Riekd l)2002|) describe a UCHII re- 
gion associated with the IRAS point sou rce, and quote 
a dista nce of 2.2 kpc ijKomne et al.lll98'9T) . and IBik et alJ 
(2005) obtain a spectroscopic distance of 1.0-1.5 kpc. 
With no a priori reason to prefer one distance over the 
other, and no uncertainties associated with either, we use 
the average distance of 1.5 kpc for our analysis. 

3.2.1. The KLF 

We present color-color and color-magnitude diagrams 
for this cluster in Figures and 03 The "cluster region" 
was defined to coincide with the extent of the near-IR 
nebulosity, and field stars were statistically corrected as 
described in 13.11 The CMD shows objects spanning a 
range of extinctions, with relatively few objects with col- 
ors consistent with unextincted main-sequence stars re- 
maining in the statistically corrected data. Since the 
cluster in this case did not fill the field of view, we were 
able to use the data from the non-cluster portions of the 
region as our field, eliminating the need for an off-source 
2MASS field and extrapolation based on the KLF. The 
J — H vs. H — K color-color diagram shows that the ma- 
jority of sources in the cluster region exhibit if -band ex- 
cess and fall in the region populated by reddened CTTS 
and YSOs, suggesting they may be pre-main-sequence 
objects. Because the amo unt of -Ff-band excess is affected 
by ma ny factors (see, e.g. lMever. Calvet. fc Hillenbrand! 
( 1997)), deriving masses for these objects is difficult. We 
thus derive only a KLF for this cluster, and do not con- 
vert it to an IMF or derive an IMF from the CMD. We 
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determine and correct for our incompleteness as in 13.11 
The KLF we derive using a statistically-corrected sam- 
ple of all sources detected in K, with no attempt to cor- 
rect for extinction due to the uncertainty of the intrinsic 
H — K color of the YSOs that are present, has a slope of 
0.30 ± 0.03. If only sources detected in H are included, 
the KLF declines for K > 14.5 since the high extinction 
means the fainter objects are less likely to be detected in 
H. 

3.2.2. Pre- Main- Sequence Objects 

A total of 37 out of 58 sources (63%) detected in all 
three bands show a K-band excess in the color-color di- 
agram, suggesting they are pre-main sequence objects. 
This is a lower limit on the actual pre-main-sequence 
fraction of the cluster, since objects with a sufficiently 
high IR excess may be detected in K but not in J or H. 
A total of 49 sources were detected in K within the clus- 
ter region that were undetected in J, H, or both. Adding 
in these sources would give a PMS fraction of 80%. The 
latter figure is an upper limit, since some of the K — only 
detections are likely to be knots of nebular Br 7 emission 
or heavily extincted background stars. Comparing these 
values to the near-IR excess frac tion of embedded clus- 
ters o f known ages presented in lHaisch. Lada. fc Ladal 
(|2001t) . we conclude that the age of the IRAS 6058 clus- 
ter is less than 3 Myr. 

We observe a few sources with colors even red- 
der than the reddened exten si on of the CTTS locus. 
iMever. Calvet. fc Hillenbrand! l)1997ft observe sources 
with similar colors, and suggest re-radiation by an ex- 
tended envelope as an explanation. 

3.3. IRAS 06104+1524 Cluster 

The near-IR cluster image (Figure|5J) of the cluster near 
IRAS 06104+1524 shows a clear separation into two sub- 
clusters separated by ~ 2'. The southwest subcluster is 
dominated by two closely spaced bright sources while the 
northeast subcluster is denser and is not dominated by a 
single object. A ridge of marginally higher density than 
the surrounding field appears to lie between the subclus- 
tcrs, though it is not apparent whether this is a real 
feature. The MSX 8 /im image similarly shows two sep- 
arate peaks, with no indication of a connection. These 
are treated as separate clusters bv lBica et alJ l)2003l) , and 
there are IRAS point sources associated with each of 
them (IRAS 06104+1524 and IRAS 06103+1523, respec- 
tively). The IRAS p oint sources both have kin ematic 
distances of 3.5 kpc ijWouterloot fe Brand! IT9891. sug- 
gesting the subclusters are related. The radio kinematic 
distances derived to the two sources are the same, and 
the angular separation is comparable to the size of either 
subclump; in addition, a slight overdensity of stars can 
be seen in the K image. These factors suggest that, at 
the very least, these clusters are related; they may differ 
in age, but they are likely to be part of the same general 
star- formation event. We see no difference between the 
two apparent in either the CMD or the color-color dia- 
gram (Figures 1101 and ITTf) : if they do differ in age, it is 
beyond the ability of our data to discern. Due to this ap- 
parent association and the small number of stars in each 
cluster, we analyze the two together as a single cluster. 
The CMD after statistical correction and adjusting to a 



distance of 3.5 kpc l)Wouterloot fc Br and 1989) is shown 
in Figure El 

Using the color-magnitude diagram based method of 
deriving an IMF as described above, for a limiting extinc- 
tion of Ay = 25, we find an IMF slope of T = -0.9+0.25. 
Using the KLF method (with no extinction correction, to 
mimic the results of a study with only single-color pho- 
tometry available) we arrive at T = —2.6 ± 0.3. Even 
after allowing for the errors to be larger than quoted due 
to uncertainty in the photometry and the conversion to 
mass, these two slopes are inconsistent with each other, 
suggesting that systematic effects in one or both methods 
dominate over the statistical errors. 

3.4. Sh 2-288 

The near-IR cluster image (Figure of the cluster 
near Sh2-288 shows a cluster with a dense core, crossed 
near the center by a dust lane. The center of the clus- 
ter is unresolved in our images, taken with a seeing of 
0.7". This region was previously iden tified as an em- 
bedded cluster by (D utra fc Bicall2001[). In thei r cata- 
log of outer-Galaxy HII regions, iRudolph et all 1(1991 
quote widely disparate distances, with a radio kinematic 
dist ance of 7.2 kpc and a photometric distance of 3.0 
kpc l|Brand fc Blitzll993(l . The kinematic distance would 
make the sources we observe (several with K < 12) ex- 
tremely massive, it is thus far more likely that the pho- 
tometric distance is correct, and we have used the pho- 
tometric distance for our analysis. 

The 8 /im image from the MSX mission shows a peak 
coincident with the near-IR nebulosity; there is not a 
significant amount of 8 /im emission from regions dark 
in the NIR. This suggests that there are not a significant 
number of sources so deeply embedded that they cannot 
be seen in K present in this cluster. 

The K versus H — K color-magnitude diagram of the 
cluster near the HII region Sh2-288 (Figure I14JI clearly 
shows the effects of variable extinction; the stars sep- 
arate into two groups, one nearly unextincted and one 
with ~ Ay — 5. We correct for field star contamination 
using the region of the field outside the cluster region as 
described above. The J — H versus H — K color-color 
diagram (Figure ^J) shows few stars separated from the 
main sequence locus by more than 2er, suggesting that 
most stars in this cluster are on the main sequence. Ex- 
treme outliers in the color-color diagram were inspected 
individually; in general, they lie in the crowded central 
region of the cluster and most likely suffer from poor 
photometry due to the different PSFs in H and K that 
resulted from variations in seeing. Such sources were 
excluded from analysis of the KLF and the IMF. Addi- 
tionally, the brightest source in the cluster, which lies 
in the most crowded central region and has a FWHM 
slightly broader than most sources in the field, was ex- 
cluded since it is quite likely to be a blend of multiple 
sources. We consider that the effects on the IMF are 
likely to be worse if a blend is included than if any single 
star, even the most massive, is excluded. 

Using the photom etric distance of 3.0 kpc from 
(Rud olph et al. 1996), we derive an IMF from the KLF of 
r = —1.95 + 0.62. To better compare the two methods of 
deriving the IMF, we included only those sources which 
were also detected in H, so that the same dataset would 
be used for both the KLF and CMD methods of deriv- 
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ing the IMF. We individually de-reddened sources in the 
H — K CMD until they were on the main sequence, im- 
posing an extinction limit as before, and derive an IMF 
of r = —1.62 ±0.65. Given the large uncertainties, these 
results are entirely consistent. The better agreement may 
be because the extinction bias is less in the latter case. 

3.5. Comparison of Methods for IMF Determination 

A summary of the IMFs derived for the three clusters 
without significant number s of pre-main-seque nce stars 
and the similar results from l|Leistra et alJ2005j) is shown 
in Table In each case, there is a significant difference 
between the IMF derived from the KLF and that derived 
from the CMD. Does this reflect only uncertainties, or is 
one method in general more reliable than the other? A 
simple analysis would suggest that the CMD method is 
more reliable, simply because it uses more information; 
the extinction clearly varies across many embedded clus- 
ters (of those analyzed here, most notably Sh 2-288 and 
IRAS 06058+2158), and accounting for this should pro- 
vide a more robust estimate of the true IMF. We cannot 
say for certain that this is the case, however, without 
obtaining spectra for most of the stars in each cluster, 
so that we can classify them spectroscopically and ob- 
tain individual masses. We note that in each case, the 
IMF we derive from the CMD by individually correcting 
the extinction for each object is flatter than that we de- 
rive from the KLF by assuming a single extinction for 
the entire cluster. This suggests that more massive stars 
may preferentially lie in more heavily extincted regions 
in embedded clusters. Resolving this seeming discrep- 
ancy would require spectra for a large number of cluster 
members, so that the IMF derived from spectral classifi- 
cation of stars can be compared to that derived via the 
different photometric methods. 

We examine the relation between derived mass and ex- 
tinction (with the extinction limit imposed) for the Sh2- 
217 cluster, where we have the most data, in Figure 
Such a relation appears to be present, albeit at low sig- 
nificance. This effect is opposite in sign to what would 
be expected from massive stars clearing their immediate 
environment more rapidly than their lower-mass coun- 
terparts, but the stars we observe are mostly of interme- 
diate mass, rather than truly high mass, such that their 
winds are not as significant; the effects of mass segrega- 
tion, placing the more massive stars at denser regions in 
the cluster, appear to dominate over the effects of clear- 
ing in these clusters. Since the clusters are numerically 
dominated by low-mass stars, the average extinction will 
be mostly determined by the average for the low-mass 
stars, changed slightly by the average for high-mass stars. 
If more massive stars are indeed preferentially found at 
higher extinction as our results suggest, this would mean 
the majority of (low-mass) stars are over-corrected for 
extinction by a small amount when using a single value, 
while a few (high-mass) stars are under-corrected (thus 
lowering the derived mass) by a large amount; thus, the 
effects of over-correcting and under-correcting extinction 
do not fully cancel out, since the large under-corrections 
would be more likely to move stars between mass bins 
than the small over-corrections. If there is no relation 
between mass and extinction, we would expect these two 
results to cancel, since the average extinction for low- 
mass stars would be the same as that for high-mass stars. 



4. SUMMARY 

We present NIR images of four embedded clusters in 
the outer Galaxy. In the case of the cluster near IRAS 
06058+2158 the number of stars with NIR excess indi- 
cates a pre-main-sequence fraction between 60% and 80% 
and an age of less than 3 Myr; the other three clusters 
show less nebular emission and fewer stars with NIR ex- 
cess indicating an older age. We compute the IMF for 
the three clusters dominated by main-sequence stars, in 
each case using both a KLF-bascd method relying on a 
single extinction value for the cluster and using only K 
band data and a CMD-based method where an individual 
extinction value is calculated for each star. We found a 
statistically significant difference between the two values 
in two of the three cases, prompting us to examine IMF 
values of embedded clusters from the literature to deter- 
mine whether systematic effects are at work. We found 
a significant difference in the mean value of the IMFs for 
embedded clusters derived from methods that handle ex- 
tinction individually compared with those that adopt a 
single value for the extinction. Although a larger sample 
would help to make this claim more r obust, since many o f 
the results come from a single study ijPorras 
and methodological details of that work could affect the 
results, we consider it to be significant enough that IMFs 
obtained by different methods should not be compared 
in an attempt to search for variations in the IMF from 
region to region. 

Truly reliable IMFs for embedded clusters will most 
likely require spectra for a large number of stars in the 
clusters; we are continuing to try to obtain spectra for 
these sources to better characterize the massive star pop- 
ulation and the IMF of these clusters. 
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TABLE 1 

IMFS FOR EMBEDDED STELLAR CLUSTERS 



Cluster 


IMF: Common Ay 


IMF: Individual A v 


Distance (kpc) 


Source 


Sh2-288 


T = —1.95 ± 0.82 


T = —1.62 ± 0.5 


3.0 


This work 


IRAS 06104+1524 


T = -2.49 ±0.3 


T = -1.38 ±0.6 


3.5 


This work 


Sh2-217 


T = -2.7 ±0.25 


T = -1.61 ±0.2 


5.0 


This work 


G305.3+0.2 


T = -1.5 ±0.3 


T = -0.98 ±0.2 


4.0 


Paper I 


- 


i 


■ 
■ 








Fig. 1. — iiT-band image of the Sh 2-217 cluster center (North=up, East=left). Image is approximately 120" on a side. The stellar 
density docs not plateau in our entire 3'FOV, suggesting the cluster outskirts continue at least to the edge of our image. 
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Fig. 2. — K vs. H — K color-magnitude diagram for the cluster near Sh2-217. Average DAOPHOT error bar is the size of the plot 
symbols or smaller. Overall uncertainties including calibration (which include terms that will not affect the relative position of the points) 
are indicated by the symbol in the upper right. 
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Fig. 3. — Statistically-corrected K vs. H — K color-magnitude diagram for the cluster near Sh2-217. Statistical correction was done 
based on the 2MASS Point Source Catalog for an adjacent field, extrapolating to the our limiting magnitude based on the luminosity 
function. 
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Fig. 4. — K luminosity function for the Sh2-217 cluster. (Solid is measured; dashed is corrected for incompleteness.) Since the cluster 
fills the field of view and the comparison 2MASS data do not go as deep, we do not show a field sample for comparison. 
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Fig. 5. — K vs. H — K color-magnitude diagram for the Sh2-217 cluster adjusted to a distance of 5 kpc. The ZAMS converted to 
observed quantities as described is overplotted. 




Fig. 6. — Color composite (J = blue, H = green, K = red) of the IRAS 06058+21 cluster region (North=up, East=left). Image is 
approximately 140"on a side. 
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Fig. 7. — Statistically field-star corrected K vs. H — K color-magnitude diagram for the cluster near IRAS 06058+2158. Significant and 
variable extinction is evident in this cluster. Typical errorbars are the size of the plot symbols. Overall uncertainties including calibration 
(which include terms that will not affect the relative position of the points) are indicated by the symbol in the upper right. 
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Fig. 8. — J — H vs. H — K color-color diagram for the cluster near IRAS 06058+2158. The lines (parallel to the reddening vector) 
delineate the possible location of reddened main-sequence stars. Over half of our sources show if -band excess; their colors are not consistent 
with reddened main-sequence sources. Typical errorbars are the size of the plot symbols. Overall uncertainties including calibration (which 
include terms that will not affect the relative position of the points) are indicated by the symbol in the upper right. Points to the left of 
the reddening lines lie in crowded regions and may suffer from confused photometry. Points more than 3 a to the left of the reddening line 
were excluded from the IMF determination. 
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Fig. 9.— Color composite (J = blue, H = green, K = red) of the IRAS 06103+1523 / 06104+1524 cluster(s) (North=up, East=left). 
Image is approximately 120"on a side. 
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Fig. 10.— K vs. H-K color-magnitude diagram for the IRAS 06104+1524 / 06103+1523 cluster(s) adjusted to a distance of 3.5 kpc. 
Circles: NE cluster. Asterisks: SW cluster. There is no difference in the CMD apparent for the two clusters, so we treat them as one to 
improve the statistics. Typical errorbars are the size of the plot symbols. Overall uncertainties including calibration (which include terms 
that will not affect the relative position of the points) are indicated by the symbol in the upper right. 
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Fig. 11— J - H vs. H-K color-color diagram for the IRAS 06104+1524 / 06103+1523 cluster(s). Typical errorbars are the size of 
the plot symbols. Overall uncertainties including calibration (which include terms that will not affect the relative position of the points) 
are indicated by the symbol in the upper right. 
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Fig. 12.— Statistically corrected K vs. H-K color-magnitude diagram for the IRAS 06104+1524 / 06103+1523 cluster(s) adjusted to 
a distance of 3.5 kpc. Typical errorbars are the size of the plot symbols. Overall uncertainties including calibration (which include terms 
that will not affect the relative position of the points) are indicated by the symbol in the upper right. 




Fig. 13. — Color composite (J 



= blue, H = green, K = red) of the Sh 2-288 cluster region (North=up, East=left). 
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Fig. 14. — Statistically corrected K vs. H — K color-magnitude diagram for the Sh2-288 cluster. The bright central source, which 
we have identified as a blend of two or more stars, has been removed as well. Typical errorbars are the size of the plot symbols. Overall 
uncertainties including calibration (which include terms that will not affect the relative position of the points) are indicated by the symbol 
in the upper right. 
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Fig. 15. — Statistically corrected J — H vs. H — K color-color diagram or the Sh2-288 cluster. The lines (parallel to the reddening 
vector) delineate the possible location of reddened main-sequence stars and of reddened T Tauri stars. Typical errorbars are the size of the 
plot symbols. Overall uncertainties including calibration (which include terms that will not affect the relative position of the points) are 
indicated by the symbol in the upper right. 



Embedded Stellar Clusters: II 21 




Fig. 16. — Mean extinction derived from H — K color for sources in each mass bin for the Sh2-217 cluster. An extinction limit of Ay = 10 
has been imposed; sources at higher derived extinction are not included in the calculation. More massive sources appear to preferentially 
lie at higher extinctions, albeit with low significance. 



